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ABSTRACT 

This  paper  presents  tables  and  graphs  assessing 
the  accuracy  of  a formula  (programmable  on  a desk 
calculator)  for  approximating  the  circular  coverage 
function  of  a circular-normal  probability  distribution. 
Use  of  this  formula  together  with  the  error  curves 
presented  here  enable  one  to  calculate  the  circular 
coverage  function  with  a maximum  error  on  the  order 
of  0.0002  for  offsets  up  to  5o  and  covering-circle 
radi i up  to  k.5o. 


ON  APPROXIMATING  THE  CIRCULAR  COVERAGE  FUNCTION 


INTRODUCTION 

Consider  a circular  normal  probability  distribution  with  standard 
deviation  a centered  at  a point  A.  The  probability  of  falling  within 
a distance  R of  a point  B r units  away  from  A is  known  as  the 
circular  coverage  function  P(R/a,  r/a) . It  is  given  analytically  by 
definite  integrals  or  infinite  series  that  cannot  be  expressed  in  closed 
form: 
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(reference  1),  where  R'  = R/o,  r'  = r/a,  and  1^  is  the  modified  Bessel 
function  of  order  n. 

EVALUATION  OF  THE  FUNCTION 

Extensive  tables  of  this  function  are  available  (references  1 through 
1*  and  other  references  cited  in  reference  5).  However,  it  is  more  convenient 
to  use  closed-form  approximations  that  can  be  programmed  on  a desk  calculator 
if  the  accuracy  desired  is  not  too  high.  A standard  reference  handbook 
(reference  6)  gives  the  following  approximation  of  the  circular  coverage 
function: 
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rf(R'),  if  R*  <1 
P(R',  r')  - { 

U - Q(x),  if  R'  > 1 

where 


a(x)  = (2Tr)'’^^  f e‘‘ 

X 


The  formula  for  P can  be  programmed  on  a TI-59  calculator  using 
127  steps  (table  1).  The  approximation  used  for  the  normal  distribution 
Q(x)  is  (reference  6): 


2 

Q(x)  = (2tt)"’^^  e"’'  (b,t  + b2  t^  + b^  t^  + b 4 t^  + b^  t^)  + E(x) , 

_8 

where  the  error  E(x)  is  smaller  than  7-5  x iO  and 


(1  + px)  ^ 
.23161*19 
.319381530 
-.356563782 


bj  - 1.7811*77937 
bjj  = -1.821255978 

bj  = 1.330271*1*29 
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ACCURACY  OF  THE  APPROXIMATION 

Reference  6 gives  no  source,  method  of  derivation,  or  statement  on 
accuracy.  The  purpose  of  this  paper  is  to  evaluate  the  error  made  in  using 
the  approximation  formula. 

The  approximation  for  R'  < 1 is  the  well-known  Carlton  formula 
(reference  1).  For  R'  < 0.i|,  it  is  accurate  to  k decimal  places.  For 
0.1*  < R'  <1,  it  is  off  by  no  more  than  .0065  (and  usually  by  no  more  than 
.001  - reference  1). 

For  the  cases  R'  ^ 1 , direct  comparison  with  the  tables  (references 
2 and  3)  reveals  that  the  maximum  error  is  .024  for  ' ^ ^*5  and 

0 < r'  ^ 5.  The  error  is  usually  smaller. 

Figures  1 and  2 are  graphs  of  the  difference  between  the  exact  and 
approximate  values  of  the  circular  coverage  function,  with  linear  inter- 
polation used  between  the  following  milestone  values;^ 

R'  - .25,  .5  , .75,  .9,  1,  1.1,  1.25,  1.5,  1.75,  2,  2.5,  3,  3.5,  4,  4.5 
r'  = 0,  .5,  I,  1.5,  2,  2.5,  3,  4,  5. 

Figure  3 is  a graph  of  the  ratio  of  the  approximate  value  of  the 
circular  coverage  function  to  its  exact  value,  with  linear  interpolation 
used  between  the  above  milestone  values.  In  table  2 we  compare  exact  and 

'since  the  formula  exhibits  a discontinuity  at  R'  - 1,  the  graphs  for 
•9  < R'  <1  are  not  accurate. 
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approximate  values  of  the  circular  coverage  function  for  these  milestones: 
vaiues,  difference,  ratio,  and  percent  increase  of  exact  vaiue  over 
approximate  value. 

As  can  be  seen  from  table  2,  whenever  the  exact  value  is  smaller 
than  the  approximate  value,  it  Is  by  no  more  than  2.k%.  Further,  whenever 
the  exact  value  is  over  11%  larger  than  the  approximate  value,  the  exact 
value  is  smaller  than  .16.  If  over  25%  larger,  P < .08.  This  indicates 
that  It  Is  for  very  small  vaiues  of  the  circular  coverage  function  that 
the  approximate  formula  has  a higher  chance  of  being  off  by  a significant 
percentage. 

In  conclusion,  the  programmable  formula  for  the  circular  coverage 
function,  together  with  figures  1 and  2,  enables  one  to  evaluate  the 
function  with  a maximum  error  on  the  order  of  .0002  for  0 ^ 8'  _5  ^-5  and 
0 < r'  <5. 
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TI-59  PROGRAM  FOR  CIRCULAR  COVERAGE  OF  A 
CIRCULAR  NORMAL  Dl STRI BUTIONa 
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The  offset  r Is  stored  in  register  09,  the  radius  R of  the  covering 
circle  In  08,  and  the  standard  deviation  a in  10.  After  r and  a 
are  stored,  keying  in  R and  pressing  R/S  will  initiate  the  program. 
Occasionally  the  display  will  flash  due  to  underflow  in  the  calculation, 

but  the  displayed  value  will  be  correct.  The  master  library  module  must 
be  inserted. 
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FIG.  2:  ERROR  MADE  IN  USING  APPROXIMATION  TO  CIRCULAR  COVERAGE  FUNCTION 
FOR  CIRCULAR  NORMAL  DISTRIBUTION  (2  < R/a 
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COMPARISON  OF  EXACT  AND  APPROXIMATE  VALUES  OF  THE 
CIRCULAR  COVERAGE  FUNCTION  OF  A CIRCULAR  NORMAL  DISTRIBUTION 


is  offset,  R is  radius  of  coverage  circle,  and  a is  standard  deviation. 
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Chicago,  llbnlos,  29  September  1977.)  Oct  1977,  6 
PP- 

PP206 

McConnell,  James  M.,  “Strategy  and  Miasioiu  of  the 
Soviet  Navy  In  the  Year  2000,"  48  pp.,  Nov  1977, 
(To  be  presented  at  e Conference  on  Probleme  of 
Sea  Power  as  we  Approach  the  21st  Century,  spon- 
sored by  the  American  Enterprise  Institute  R>r 
Public  Policy  Research.  6 October  1977,  and  tub- 
aequentty  published  fat  a collection  of  papen  by  the 
Institute) 
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Pr207 

Goldbtrg.  Uwfenc*.  **Coit<Erfectivcnen  of  Po- 
tmtlal  Federal  Pobciee  Affeccuig  Reeeerch  A 
Dewlopment  Expendituree  In  the  Auto.  Steel  and 
Food  Induatrkt,^  36  pp..  Oct  1977.  (Presented  at 
Southern  Economic  Association  Meetings  beginning 
2 November  1977) 

pp2oe 

Roberts.  Stephen  S..  '*The  Decline  of  the  Overseas 
Slatiofi  Fleea:  The  United  States  Aiiitie  Fleet  and 
tte  Shanghai  CrMs,  1932.**  18  pp..  Nov  1977.  (Re. 
printed  from  The  American  Neptune,  Voi. 
XXXVII..  No.  3.  July  1977) 

PP  209  - Classified. 

PP210 

Kassing,  David,  "Protecting  The  Fleet,**  40  pp..  Dec 
1977  (Prepareo  r the  American  Enterprise  Insti- 
tute Confereirce  on  Problenu  of  Sea  Power  as  We 
Approach  the  21st  Century,  October  6-7,  1977) 

FP  211 

Mizrahi,  Maurice  M..  "On  Approximating  the  Cir- 
cular Coverage  Function.**  14  pp.,  Feb  1978 

PP2I3 

Mangel.  Marc.  "Fluctuations  in  Systems  with 
Multiple  Steady  States.  Application  to  Lanchesler 
Equations,**  12  pp..  Feb  78.  (Presented  st  the  First 
Annual  Workshop  on  the  Information  Linkage  Be- 
tween Applied  Mathematics  and  Industry,  Naval  PC 
School.  Feb  23-25.  1978) 

PP216 

Coile.  Russell  C..  "Lotka’i  Frequency  Distribution 
of  Scientific  Productivity.**  18  pp..  Feb  1978.  (To 
be  tubmiffed  for  publiceiion  m American  Society 
for  Infomution) 
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